Is there anyway in which this algorithm can be expressed in terms of canonical matrix operations.
2.
Where E ^ { \ gamma } _ { \ delta } ( \ alpha \ beta ) is the canonical matrix with all components equal to zero except for the entry corresponding to row \ alpha and column \ beta, that is equal to one.